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Abstract: We investigate supersymmetric models for dark matter which is represented by 
pseudomoduli in weakly coupled hidden sectors. We propose a scheme to add a dark matter 
sector to quiver gauge theories with metastable supersymmetry breaking. We discuss the 
embedding of such scheme in string theory and we describe the dark matter sector in terms 
of D7 flavour branes. We explore the phenomenology in various regions of the parameters. 
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1. Introduction 

Cosmological observations have established the existence of dark matter which is not com- 
posed by any of the Standard Model (SM) particles, and with a relic abundance of the 
order of Vlh? ~ 0.1 (for reviews see Jl|, ||] and reference therein). 

A stable particle in thermal equilibrium with the SM in the early universe is usually 
referred as cold dark matter (DM). As long as the universe expands, the DM ceases to 
annihilate efficiently and freezes out, leaving a relic abundance Q 

ntf „ _L „ E**l (i.i) 
M m 2 dm 

where h is the Hubble parameter in units of 100 km/sec per Mpc, (av) is the annihilation 
cross section, a/jM is the characteristic coupling and M^m is the mass of the dark matter 
particle. For a massive particle with weak interaction {goM ~ 1); i-e. a WIMP, the required 
relic abundance is obtained for Mdm ~ 1 TeV. The TeV scale that naturally appears 
suggests some relations between DM and electroweak symmetry breaking (EWSB), and 
has been dubbed as the WIMP miracle. 

A standard explanation for the origin of the electroweak scale is supersymmetry break- 
ing (for a phenomenological introduction to supersymmetry and references see ||). The 
idea that supersymmetry and its breaking should account also for the origin of the DM 
has prompted an extensive investigation. Much effort has been focused on the mechanism 
of gravity mediation, where the lightest neutralino of the MSSM should be a stable, viable 
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DM candidate. A phenomenological drawback for models of gravity mediation is that they 
do not address the susy flavour problem. 

The flavour problem is solved in the gauge mediation scenario [0, |], because the SM 
gauge interactions are flavour blind. In gauge mediated models supersymmetry is broken 
at low energies and the gravitino is the LSP. However, the gravitino is typically too light 
(with mass in the range of eV to Gev) |], 10 to be a viable cold DM candidate. Other 
alternatives for susy DM have been and are under inspection (see for instance |11, |l^, |lS|]). 

In gauge mediated theories particles in the hidden sectors can be alternative DM can- 
didates. In such models, supersymmetry breaking in the hidden sector occurs dynamically, 
via strong dynamics effects. The DM can then be identified with the mesons and the 
baryons of the strongly coupled hidden sector [ 14 1 . 

On the other hand, one can explore the possibility of realizing dark matter in weakly 
coupled hidden sector with spontaneous supersymmetry breaking IS, 16]. These models 
are relevant since they can arise as the low energy description of UV free gauge theories, 
as shown by [jlTj . 

Such weakly coupled models typically break supersymmetry in metastable vacua, with 
a spectrum of elementary particles. The scalar potential often presents tree level flat 
directions, which are lifted by one loop quantum corrections. These pseudomoduli fields 
have weak scale interactions and their one loop mass is of the order of the supersymmetry 
breaking scale, i.e the TeV scale naturally arises. If they are stable against decay because of 
a discrete Z2 symmetry, they represent viable cold dark matter candidates. This possibility 
has been explored in O'Raifeartaigh like models in [|l^] and also analyzed in pq| . 

In this paper we shall study models of pseudomoduli as dark matter and wish to answer 
some questions concerning the relative UV completion. This is a non trivial problem, with 
severe constraints [18, [0| for consistent embedding in a UV complete theory. We shall 
propose a quiver model for which a stringy origin can be reached. We embed models with 
pseudomoduli DM in appropriate quiver gauge theories which arise from D-branes at CY 
singularities. The starting point is a quiver theory, inherited and/or interpreted as an IR 
Seiberg dual, which provides the hidden sector of dynamical breaking of susy. The next 
step is the addition of a dark matter sector, represented by an extra node in the quiver 
connected through matter interactions to the hidden sector. The addition of D7 flavour 
branes in the CY can be put in correspondence with the dark matter sector in the quiver. 

The structure of the paper is the following. In section |2| we propose the general strat- 
egy to embed pseudomoduli dark matter in quiver gauge theories, we comment on the 
phenomenological constraints and also on the string theory realization. In section ||| we 
build a concrete example; we couple pseudomoduli DM to the KOO model [20], and we 
discuss the related phenomenology. In section || we provide for a UV origin to the model 
in terms of a step of Seiberg duality: the UV model is obtained by deforming an L 131 non 
isolated singularity and by adding flavour D7 branes. A conclusion follows. In appendix [A| 
we review a basic cosmological bound on the supersymmetry breaking scale. In appendix 
[B] we review the procedure of flavoring with D7 branes and we discuss its relation with 
the DM sector. In appendix |C] we provide the details of the one loop computations for the 
pseudomoduli masses. 
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Figure 1: Quiver representing the dark matter coupling to the supersymmetry breaking sector 



As we were finishing this paper, we were informed of pT] which study leptophilic dark 



matter [22] in quiver gauge theories. 



2. Pseudomoduli DM in quiver gauge theories 



As anticipated, quiver gauge theories are useful settings for the study of the stringy origin 
of pseudomoduli dark matter. In the quiver, we distinguish between a supersymmetry 
breaking sector and a DM sector. The first one is thought as an ISS like model: it is 
characterized by a metastable vacuum, a weakly coupled SU (N) gauge symmetry and an 
SU(F) flavour symmetry. We also require i?-symmetry to be broken in this sector. We 
parametrize this sector with a spurionic chiral field X, an SU(N) singlet, which acquires 
a vev and an .F-component, i.e. X = M + 6 2 F, and its fermionic component is the 
Goldstino. The chiral field X couples to two chiral fields Q12 and Q21, respectively in the 
fundamental and antifundamental representation of the gauge group SU(N), and in the 
antifundamental and fundamental of the flavour symmetry SU(F). A three-linear coupling 
XQ12Q21 induces a mass splitting between the bosonic and fermionic components of the 
fields Q12 and Q21, which are the messengers of susy breaking. This sector is depicted in 
the rectangular region in the figure [l|. 

We then add the DM sector, characterized by a U(l)d gauge symmetry and two pairs 
of bifundamental fields (see figure |l|). These fields interact with the messengers through 
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the superpotential 

Wdm = QviQ2dQd\ + Qi\Q\dQd2 + r md2Q2dQd2 (2-1) 

where the traces are understood. This new sector does not change the solution of the 
equations of motion for the hidden sector, and the metastable vacuum remains a tree level 
minimum of the scalar potential. 

In the DM sector, the fields Q2d and Qd2 have a tree level mass and are stabilized at 
zero vev. They are called link fields. 

The other chiral multiplets Q\d and Qdi are massless at tree level. In particular, their 
scalar components are tree level flat directions not associated with any broken global sym- 
metry, i.e. pseudomoduli. Both the fermionic and bosonic components of these multiplets 
Qid and Qdi can get a mass at one loop. Since supersymmetry is broken, we expect scalars 



and fermions to get different masses. Typically, as discussed in [15, 16], the scalar masses 
are higher than the fermionic masses. This implies that in our setting the fermions ipQ ld 
and ipQ dl have lower masses than their scalar partners. These fermions are the cold DM 
candidates. 

Brief and comments 

The previous scheme can generate viable dark matter candidates provided the hidden sector 
fulfills some basic requirements. 

The mechanism that we consider for communicating the supersymmetry breaking to 
the visible sector (MSSM) is gauge mediation [|], H|. This is a standard technique in susy 
breaking quiver gauge theories [23, ^4|, 25]. Here we shall concentrate on the direct gauge 



mediation scheme. 

Direct gauge mediation is realized by embedding the MSSM gauge group in some 
subgroup of the hidden sector. Different choices lead to charged or to uncharged dark 
matter under the MSSM gauge group. Dark matter is charged if we identify the MSSM 
gauge group with SU(F), whereas it is uncharged if we embed the MSSM gauge group in 
SU (N). Alternatively, we can give mass to the other pair of chiral fields adding midQidQdi 
to the superpotential (|2.l|), setting rri2d = 0. This exchanges the role of dark matter and 
link fields as well as that one of charged and uncharged dark matter. In this paper we 
discuss the case of uncharged dark matter, which is less constrained by experiments. In 
this case the massive link fields are charged under the MSSM gauge group and the U(l)d 
dark matter gauge group. The dark matter is charged under the U(l)d gauge group and a 
hidden sector flavour group. The interaction between the MSSM gauge group and the dark 
matter is obtained via kinetic mixing Recently this mechanism has been investigated 

0.0- 



m 



In models of gauge mediation i?-symmetry has to be broken for the gaugino to acquire 



non trivial mass. The models analyzed in [15] were required to respect a spontaneously 



broken i?-symmetry that prevents tree level mass terms and extra couplings for the pseu- 
domoduli dark matter candidates. We also admit the explicit breaking of i?-symmetry in 
the hidden sector. Indeed in our case the UV stringy origin uniquely settles the structure 
of the interaction. 
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There should be a discrete Z2 symmetry that forbids dark matter to decay. In our 
case this symmetry follows automatically from the non-chiral structure of the DM sector 
and from its interaction superpotential fl2.ip . Furthermore the pseudomoduli dark matter 
Q2d and Qd2 have to be stabilized at one loop at the origin of the pseudomoduli space such 
that this Z2 symmetry is unbroken. 

As discussed in [|l5|, 16], we expect the one loop masses of the DM fermions to be 



smaller than the masses of their scalar superpartners. This high difference between scalar 
and fermion masses implies that the decay of the Q2d an d Qd2 scalars does not affect the 
dark matter relic density. This property of the one loop masses for the scalars and the 
fermions has to be checked in any model of pseudomoduli DM. 

We require to have a TeV scale dark matter mass, and also a not too heavy superpartner 
spectrum in th MSSM. Hence we demand the parameter 

R ee Mm (2.2) 
m\ 

to be of order 1, where we estimate the soft mass scale with respect to the gaugino mass m\. 
The dominant contribution to the dark matter cross section comes from its annihilation 
into dark photons 

M " TFT- ( 2 - 3 ) 

1V1 DM 

where ad is the coupling of the U(l)d gauge group and Mdm ls the dark matter mass. For 
gd ~ 1 and Mdm ~ O(TeV) this cross section secures a satisfactory relic abundance, and 
it avoids dark matter overabundance |15L Qq| . Note that for this annihilation to be efficient 



the dark gauge boson mass needs to be lower than the dark matter mass [27, [0|]. We will 
not address the naturalness of this new scale here. 

The realization we provided of pseudomoduli DM in gauge theories is rather generic. 
Many models of supersymmetry breaking in metastable vacua with explicitly broken R- 
symmetry have been realized [28, |30|, |3^, 33, 34, |35], 36]. The procedure for adding 
a DM sector just explained can be applied to all these models. 

The stringy origin of the DM sector 

Metastable supersymmetry breaking is common to gauge theories arising from D3-branes 
at CY singularities. In this framework the DM sector corresponds to .D7-flavour branes. 
Indeed, by properly adding Z?7-branes, we obtain the field content and the interaction 
superpotential of the DM sector. In the appendix ^ we discuss the realization of DM sector 
as D7-branes. This suggests a connection between pseudomoduli DM and D7-branes at 
CY singularities. 



3. Coupling DM to the KOO model 

In this section we give a concrete example of the strategy presented above. We consider 
as the supersymmetry breaking sector the KOO model I^Ofl . The low energy description of 
this theory is a three node quiver gauge theory U(N{) x U{N2) x U{N%), where 

N 1 = N 2 = N N 3 = N + M (3.1) 
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Figure 2: Quiver representing the flavored KOO model 



We choose N < M such that the [/(iVjj) gauge group is infrared free. The other gauge 
groups are considered as very weakly coupled at low energy. The superpotential for this 
model is 

Wroo = -hulXu - hf4X 33 + hm 13 X n X 3 3 + hq-X -q (3.2) 



where 




V A 31 A 33 / 



This model breaks supersymmetry in a long living metastable vacuum and does not have an 
i?-symmetry. The absence of .R-symmetry implies that direct gauge mediation is a viable 
mechanism to transmit supersymmetry breaking to the superpartners in the MSSM ]^(| . 

We now add the dark matter sector. As already explained we add a new gauge group 
U(l)d with bifundamental fields connected with the groups U{N2) and U(N 3 ). The result- 
ing quiver is depicted in figure ||. The new fields interact with the supersymmetry breaking 
sector fields via the superpotential 

W DM = hq 23 q 3d q d2 + hq 32 q 2d q d3 + hm 2 q 2d Qd2 (3.4) 

The metastable supersymmetry breaking vacuum of the KOO model is not destabilized by 
the deformation ( |3.4j ). Hence the vacuum expectation values of the fields are 

(3.5) 
(3-6) 

The vev of the fields q± 2 and q 2 i break the groups U(Ni) x U (N 2 ) to the diagonal subgroup 
U(N)i_ 2 . The fields \i Y, Y are pseudomoduli. Their scalar components get masses via one 
loop corrections. The Goldstino is a mixture of the fermionic component of x an d of X\\ 




|3q , |37f , and it is eaten by the gravitino in the super Higgs mechanism [38]. The gravitino 
is typically too light to be a viable dark matter candidate in this gauge mediation scenario. 
The fermionic components of the fields Y, Y are then natural dark matter candidates. We 
shall perform a quantitative analysis of the spectrum to justify this scenario. 
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As a standard procedure we expand the superpotential around the vacuum 

<?3d = y 9d3 = Y Q2d = $5 Qd2 = $6 (3.8) 
and then study the quantum infrared model which is a generalized O'Raifeartaigh model 

W = hx^i^2 ~ hulx + fyUl($l$4 + Ml*2$3) + /imi 3 $3$4 

+ /iy$i$ 6 + /lY$ 2 $5 + /im 2 $5$6 (3.9) 



The pseudomodulus x get one loop corrections only from the first line in (3^9), which is 



the same microscopic superpotential as in [20]. The second line in (3J)) is the dark matter 



sector. The massive fields <3?5 and $6 are the link fields. The scalars and the fermions of 
the fields Y and Y get both one loop masses. In the appendix |C] we give the analytical 
calculations of these masses. Here we discuss the phenomenology of the model. 

Phenomenology 

The stability of the metastable supersymmetry breaking vacuum requires 

Ml > M3 Mi > rni3 (3.10) 

where fii is the messenger mass scale, /U3 is the supersymmetry breaking scale, and mi3 
is the /^-symmetry breaking mass. As already explained, supersymmetry breaking is com- 
municated to the MSSM via direct gauge mediation. The GUT SU (5) gauge group can be 
embedded both in U(N^) or in the diagonal subgroup U(N)i-2- In the first case the pseu- 
domoduli DM is charged under the GUT group. We will not investigate this possibility. 
We choose the second case, leading to uncharged DM. The gaugino mass is ]20| 

16vr 2 fij V Mi 

The scalar masses are of the same order provided that mis ~ Mi/V-^V" + M. 

We now discuss the quantum aspects of the dark matter sector. The scalar component 
of the chiral fields Y, Y acquire positive squared masses and are stabilized at the origin 
of the moduli space. This preserves the Z2 discrete symmetry that makes the DM stable. 
The DM, i.e. the fermionic component of the field Y, Y, also acquire one loop masses as 
well. 

A detailed computation of the 1-loop scalar and fermionic masses for the components 
of the chiral fields Y and Y , at all order in the supersymmetry breaking scale, is carried out 
in the appendix Here we report the analytic result at the combined third order in the 
adimensional parameters m.13/ Mi, W2/M1 an d x/Ml> an d at first order in the supersymmetry 
breaking scale The fermion mass results 

= { N + M)-^- 2 ^(l + l^f- 1 -4) + --- (3-12) 
Y 16tt z fif 3 /if 3 /j,f 



-a = ^(N + M)^H +0 r|] (3.11) 
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Figure 3: Ratio m^ Y ^ 9 /my where with my we denote the lowest eigenvalue of the scalar mass 
matrix. The masses are evaluated at the first order in the susy breaking scale, such that it cancels 
out in the ratio. We fix the couplings to h = 1, g = 1. 



The scalar diagonal and off diagonal masses are respectively 



m YY * = m Yy , = (AT + M)— -^(1 - -4- - 2^) + . . . (3.13) 



"4f = -^ + ^^ + - (3-14) 

The eigenvalues of the scalar masses can be obtained diagonalizing the resulting mass 
matrix. The off diagonal components are subleading and hence the main contribution to 
the eigenvalues comes from the diagonal masses ( 1H3|) . In the expressions above we should 
insert the vev of \ as a function of the other parameters, which is found by minimizing the 
effective potential. 

Note that there should be at least one order of magnitude between the lowest eigenvalue 
of the scalar mass matrix and the fermion mass, otherwise the DM relic abundance is 
affected by the decay of the scalars Y and Y. In the appendix |A] we review how this 
constraint gives a bound on the scale of susy breaking. In figure |3] we plot the ratio of 
fermion and the lowest scalar mass as a function of the parameters of the model: it shows 
that there are two order of magnitude between the two masses. This constraints the susy 
breaking scale to /U3 < 10 5 TeV (see appendix |A|). 

As already explained, we require the DM to be of the same order of the soft masses 
(~ 1 TeV). We estimate the parameter R fl2.2| ) as 



R=^^^^ + ... (3.15) 
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Figure 4: Ratio R as a function of the parameters TO13//Z1 and 77124/ Vi- 



and we plot it in figure ||. Figure |] shows that we can find a range in the parameter space 
where the ratio R is of order 1. Notice that the enhancement factor (N + M) cancels in 
the ratio R. This is a specific feature of the model, since the messenger fields and the dark 
matter are charged under the same global symmetry U(N^). 

As usual in theories with direct gauge mediation, this model can suffer from a Landau 
pole problem. As pointed out in p7| , requiring perturbative unification of the couplings 
below the Landau pole forces the messenger scale and the supersymmetry breaking scale to 
be large. This results in a too large mass for the gravitino, outside the cosmological bound 
worked out in [39]. We leave a detailed analysis of the issue of gauge coupling unification 



for future studies. One can solve this problem by looking at a different UV completion for 
the KOO model, for example via a cascading gauge theory. In the next section we show 
how the same low energy theory arises from a system of D3 and D7 branes probing a CY 
singularity through a Seiberg duality. 

Alternatively, one can choose a different embedding for the SM gauge group into the 
flavour group of the supersymmetry breaking sector. This embedding of SU(5)gut into 
U (-/V3) has been investigated in [37] , and it has been shown to be compatible with a gravitino 



mass which is consistent with the cosmological bound of [39|. For this embedding the 
pseudomoduli Y and Y are charged under the standard model gauge group. If we want 
to realize DM uncharged under the SM gauge group, we have to exchange the role of the 
dark matter and of the link fields. This is done by setting to zero the mass term in (|3.4|), 
and adding a new one hm^q^Qdz- The phenomenology of the low energy theory and the 
ratio R are unchanged. 
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4. Pseudomoduli DM from D-brane at CY singularities 



String theory provides a natural embedding for pseudomoduli DM. In fact, deformations of 
non isolated singularities lead to metastable supersymmetry breaking in the gauge theory 
living on _D3-branes probing the singularities |^0|, |4l| . It has also been shown |^IJ that the 
addition of flavour D7 branes can break supersymmetry in metastable vacua. Here we use 
both these effects to build a quiver gauge theory. We show that the resulting model is a 
concrete realization of the setup for pseudomoduli DM in quiver gauge theories proposed 
in section |2|. 

We consider the quiver gauge theory arising from D3 branes at the L 131 singularity. 
It can be described by the curve in C 4 

xy 3 = wz (4.1) 

and we deform it as 

x{y + Ov(y - = wz (4.2) 

We choose the rank of the four gauge group as 

AT = Ni = 1 N 2 = M N 3 = N (4.3) 
and the resulting superpotential is 

W = X11Q21Q12 + AQ12Q23Q32Q2I + m23<323<232 (4.4) 



where the mass is related to the parameter in the geometry (4.2) as m23 = £A. We add 
N4 = M — I Z)7-flavour branes associated with the fields Q 23 and Q 32 (see appendix |B|). 
This introduces new fields that interacts via the superpotential 

W D7 = PQ23Q34Q42 + PQ32Q24Q43 + rn u Q 3i Q4 3 + m 2 iQ2iQ42 (4.5) 

In figure [| we give a pictorial quiver representation of the model, with complete superpo- 
tential 

W = W + W D7 (4.6) 



We work in the range M > 2N. In this window the node 2 has the strongest coupling 
and it is in the magnetic free window. The low energy description can be obtained via 
Seiberg duality. The dual low energy theory has the superpotential 



W = Tr ( q 2 i q 2 3 <724 



/ Mn M13 M14 
M31 M 33 M 34 
^M 4 i M 43 M 44 

+ +WQ34Q34Q43 + m MQU M u Q 43 + m M Q34Q 34 M43 - /i|M 33 
where the new scales are 




+ m n M n Xn + mi3Mi 3 M 31 



(4.7) 



mn 



Ao 



"113 



AA|, 



"1Q34 = "2-34, m M QU = PA 2 , 



f4 



m 23 A 2 , nl = m 2A A 2 
(4.8) 



- 10 - 



Nl N2 N3 




N4 

Figure 5: Quiver representing the flavored L 131 model. 

Nl N2 N3 




N4 

Figure 6: Quiver representing the flavored L 131 model after Seiberg duality on node 2. 

We work in the regime where ran and mg 34 are larger than the other scales of the theory. 
We can then integrate out the fields X\i,Mu,Qs4 and Q43, and obtain the superpotential 

W = <?2lMi3<J 3 2 + g2lMi4<742 + 923^31912 + 924^41912 + 923-^33932 + 923^34942 

+ <724M 43 < ? 32 + q2AM M q 42 + m 13 M 13 M 31 + m 34 M 34 M 43 - /i§M 33 - /^M 44 (4.9) 

where we define m 34 = ^v\ [ QzJ' m Q^- Figure || is a pictorial quiver representation of the 
magnetic model. The ranks of the groups in the magnetic theory are 

Ni = 1, N 2 = N 3 = N, N A = M-l (4.10) 

We look for the metastable vacuum states. We find 

M 13 = M 31 = 0, M 34 = M 4 T 3 = 0, M U = Y, M ix = Y (4.11) 

M 33 = 0, M 44 = X, 912 = 921 = 0, 923 = 932 = M3, 924 = 942 = 

The one loop correction for the pseudomoduli are calculated after expanding the fields of 
the theory around their expectation value. Nevertheless some of their fluctuations give a 
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supersymmetric contribution at one loop. The only relevant expansions are 



942 = 4>1 924 = 4>2 912 = <p5 921 = 4>6 M 13 = (p 7 M 3 % = (p 8 (4.12) 

M 14 = Y M 41 = Y M 34 = 4 M 43 = </»3 M U = X 

Expanding around the vacuum we find the following structure for the effective superpoten- 
tial involving the pseudomoduli \-> Y an d Y 

W eff = W( X ) + W(Y,Y) (4.13) 

where 

W(x) = X4>l4>2 ~ l4x + M3(0104 + Hfo) + m-34^3^4 (4.14) 
and the superpotential for the other pseudomoduli is 

W(Y, Y) = Yfafo + Yfofc + /u 3 (<^7 + 050s) + mi 3 ^ 7 08 (4.15) 



This superpotential reduces to (3J3) in the limit 77113 > /x 3 . Indeed in this limit we can 
integrate out supersymmetrically the fields <j)-j (j>s and obtain an effective mass term for 
the fields jUg/mia </>5</>6- Hence in this limit the phenomenology of the model is the same 
as in section ||. However, also the case with mi 3 < ^ 3 is phenomenologically viable, with 
the parameter R of order 1. These features are manifest in figure [7| where we plot the 
parameter R as a function of the ratios m-i 3 //i 3 and ?n 34 //U 3 . In appendix |C] we perform 
the explicit computation for the 1-loop DM fermion mass. 



5. Conclusions 



In this article we proposed a scheme to realize pseudomoduli dark matter in quiver gauge 
theories with weakly coupled supersymmetry breaking a' la ISS. In section [2] we distin- 
guished a metastable supersymmetry breaking sector and a dark matter sector. The for- 
mer communicates the breaking of supersymmetry to the MSSM trough gauge interactions. 
This mechanism requires i?-symmetry to be broken, spontaneously or explicitly. The dark 
matter sector contains two pseudomoduli fields and it is coupled to the supersymmetry 
breaking sector through superpotential terms. This coupling induces, at one loop, super- 
symmetry breaking masses for the pseudomoduli and their fermionic partners in the dark 
matter sector. This provides the mechanism to generate a TeV scale mass for the fermions 
which are the cold dark matter candidates. 

We considered only the case of uncharged dark matter with respect to the gauge 
interactions of the standard model. There is a coupling of the U(l)d gauge group, under 
which DM is charged, to the U(1)y of the standard model via the kinetic mixing. Cases 
in which the dark matter is charged under a non abelian gauge group are left for future 
works. 

We showed that when the supersymmetry breaking sector is described by D3 branes 
probing CY singularity, the dark matter sector corresponds to the addition of D7 flavour 
branes. 
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Figure 7: Ratio R. 



In section ||| we realized the proposed scenario by coupling a pseudomoduli DM sector 
to the KOO model, and ITeV mass dark matter is rather natural. In section ||] we studied 
a model with D3 and D7 branes wrapped over a deformed CY L 131 singularity. We have 
argued that this model reduces, in the infrared, to the KOO model plus DM. 

Many extensions of our proposal can be studied. One can find other models by coupling 
the DM sector to gauge theories different from KOO. For example one can couple DM to 

, Sal M |3l, M. 



models with i?-symmetry breaking in metastable vacua [28, 29, |3(], 31 



Another interesting issue is the investigation at large of the UV completion. The generation 
of the DM sector inside a quiver gauge theory can be a result of a step of Seiberg duality, 
as we showed in the deformed and flavored L 131 theory. 
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A. Cosmological bounds 

Pseudomoduli DM arise from a tree level massless chiral multiplet Y = (4>y,^y)- It 
includes a complex scalar and a fermion. Both fields acquire masses at one loop. In the 
model we studied there is a hierarchy among these masses. The scalar mass is typically one 
or two order of magnitude larger than the fermion mass. There is a Z2 discrete symmetry 
under which Y is charged. This prevents the field ifty to decay. If it has weakly coupled 
interactions and mass at the TeV scale, it is a viable DM candidate. 

The scalar component, which is heavier, can decay in a gravitino and a fermion, <py —> 
Gipy- This decay can modify the relic density of the DM candidate ipy. We thus require 
that the decay temperature of the scalar (jjy is larger than the freeze out temperature of 
the fermionic DM ipy. The scalar decay rate is Q 

T{(Py - ^G) = <_ ( 1 _ ~ ^k- (A.l) 

The temperature associated to this decay is T ~ T 1 / 2 

Ttr _ (~y 3TeV (A . 2) 

0y VIOOTeV/ V V? / 

whereas the freeze out temperature of the fermion is Tf reeze ~ mw, y /20. The requirement 
T^ Y > Tf reeze translates in a bound on the supersymmetry breaking scale 

B. Flavoring with D7 branes 

In this paper we have proposed how to realize a dark matter sector in quiver gauge theories. 



Here we show how this new sector can be described as D7- flavour branes [42] 



We briefly review the technique introduced in |41] to add D7 branes to toric quiver 
gauge theories and to extract the interaction superpotential; we refer to the original paper 
for a detailed explanation. 

Consider a toric quiver gauge theory realized as D3 branes probing a toric CY singu- 
larity. The system can be described in terms of a dimer diagram and a useful tool is the 
Riemann surface in the mirror configuration [43]. 

By the use of these tools, it turns out that in the quiver we can associate to every 
bifundamental field a supersymmetric four cycle, which passes through the singular point, 
on which a D7-brane can be wrapped. Call one of this cycle and label the two gauge 
groups under which the bifundamental is charged as U(Ni) and U(Nj). When adding 
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Figure 8: Quiver obtained by adding a D7 brane associated with the held Xy. This is only a 
sector of an anomaly free theory. 




Figure 9: Quiver obtained adding D7 brane associated to the field Xij and X^. 

D7-branes, one should control the cancellation of RR tadpoles, corresponding, on the field 
theory side, to an anomaly free theory. 

The addition of the D7-brane adds new bifundamental fields corresponding to strings 
stretched between the D7 brane and the D3 branes. These new degrees of freedom are 
charged under the gauge groups U (iVj) or U(Nj) and under the U(1)a symmetry introduced 
by the D7 brane . We show in the figure |8| the resulting quiver. 

There is an interaction superpotential term of the type 33—37—73 that can be obtained 
analyzing the disk on the mirror Riemann surface, and it is 

W int = XijqjAqAi (B.l) 

If there are D7 branes on different four cycles, each of them introduces a couple of bifun- 
damental fields and interaction terms as in ( |B.1| ). There are also interaction terms of the 
type 37a — 7^7^ — 7^3. These terms can lead to masses for the 37 fields if a 7^7^ field 
get a non trivial vacuum expectation values. The vev breaks the U(1)a X U(1)b groups 
associated to the D7a arid D7b branes to the diagonal subgroup. 

For instance consider the addition of D7-branes associated to two different bifunda- 
mentals X^ and Xji charged under the group U(Ni) and U(Nj). The resulting theory is 

1 If there are K D7 branes on the same cycle and with the same Chan-Paton structure the symmetry is 
U(K) A . 



the quiver depicted in figure || . A vev for the 1 aIb held give raise to the following mass 
terms in the superpotential, that involves both set of flavours, 



W mass = miqAiqiB + rn 2 qBjqjA (B.2) 

in addition to the interaction superpotential 

Wi n t = XijqjAqAi + XjiqiBqBj (B-3) 

The mass parameters can be related to geometrical quantities as follows. The mass term 
corresponds in the geometry to the recombination of two D7a and D7b branes. The two 
cycles recombine in one cycle which passes at some distance e from the singular point. 
Taking local holomorphic coordinates we can parametrize the two four-cycle as z\ = and 
Z2 = 0. The mass term corresponds to a recombination of the D7a and D7b brane such 
that they now wrap the cycle z±Z2 = e. The parameter e is the distance of the four-cycle 
from the singular point, and it is related to the gauge theory parameters as 

e ~ 77117712 (B.4) 

Then the two mass parameters have to be both turned on. 

In the paper we used this technique to add flavours to quiver gauge theories and 
to extract their interaction superpotential. In section [2] we add D7-branes associated to 
the fields Q12 and Qi2- We introduced only one of the two mass terms of QB.2j ). This 
approximation can be considered as a limiting case where there is a large hierarchy among 
the two masses, and hence we neglect one of the two. This is generally the approximation 
to adopt when adding a D7-brane dark matter sector to a weakly coupled supersymmetry 
breaking sector, as in section |2] and in section ||. 

In section |] we analyzed a UV complete theory, and there we introduced _D7-flavour 



branes associated with the fields Q23 and Q32 both, with the mass terms of (B.2). The 
setup of section ^ is then obtained as the low energy description of the model by performing 
Seiberg duality. 



C. An explicit calculation 

In this appendix we show a detailed calculation of the one loop masses for pseudomoduli 
DM. We consider the superpotential 

Wx = fX + X4>i(j)2 + M<M3 + Ma) + nmfofa (C.l) 

This superpotential represents the supersymmetry breaking sector. Here i?-symmetry is 
explicitly broken by the term mi (^3 ^4. In the non supersymmetric minimum the fields <j>i 
have zero vev. The fields X is a pseudomodulus. The pseudomoduli space is tachyon free 
and stable if 

\fi 2 ± mi X\ 2 - f{ml + fi 2 ) >0 (C.2) 
The one loop analysis shows that this pseudomodulus is stabilized at (X) / |37[| . 
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DM and KOO model 

The model of section |j| is recovered by adding to ( |C.lD the superpotential for the dark 
matter sector. This is 

W 2 = m 2 5 6 + Yfafc + Yfafa (C.3) 

The fields ^5 and (f>Q are stabilized at zero vev in the non super symmetric vacuum, while 
the fields Y and Y are pseudomoduli. 

They are stabilized at one loop at the origin of the moduli space, and their scalar 
components and their fermions components get both a non zero mass. The fermions ipy 
and tpy S e t one 1°°P masses, differently from tftx, because they are not associated with the 
goldstino. However their masses are generically different from the masses of their bosonic 
partners, because they feel the effects of supersymmetry breaking. 

In this appendix we explicitly calculate the scalar and fermion masses. The calculations 



are performed by using the approach of 44]. It consists of calculating the mass term for the 
model with / turned on and repeat the same calculation, but with / = 0. The two models, 
the one with / 7^ and the one with / = 0, have the same interactions, the same field 
content, but a different spectrum. Since the masses of the fields get corrected only in the 
non supersymmetric case, the difference between the mass in the non supersymmetric case 
and the mass in the supersymmetric case coincide with the mass in the non supersymmetric 
case. Trivially we can write the equation 

m %o = m %o ~ m /=o ( c - 4 ) 

This trick reduces the number of diagrams necessary to calculate the masses of the pseu- 
domoduli. Indeed the only diagrams that contribute to the mass are the ones depending 
on / in the non-supersymmetric case. 

We first calculate the fermion mass term in the effective Lagrangian of the form 

C eff D M^ipyipY + h.c. (C.5) 
This term arises from the one loop diagrams due to the interactions 

C D 1pY(lp5<l>l + feipi) + "ipY(ip6<h + ^6^2) + h.c. (C.6) 

The calculation is not immediate, since the mass matrix for the fields (f>i and 4>2 are not 
diagonal. We have to diagonalize the bosonic mass matrices for I/O and Y = Y = 0. For 
simplicity we first diagonalize the fermionic squared mass matrix, and then we diagonalize 
the bosonic one, breaking the holomorphic structure. The eigenvalues of the fermionic 
mass matrix for the fields (f>\ , . . . , 64 are 



,2± 



ml + X 2 + 2fi 2 ± yj(m{ - \X\ 2 ) 2 + A^ 2 (m{ + 2 Re(X)mi + \X\ 



2 

The diagonal combination of the fields in the superpotential are 



(C.7) 



61= COST/91 + Sin T/?4, (f>2 = COS Tp2 Sin T /93 

63 = — sin rp2 + cos T/93 , 1^4 = — sin rp\ + cos T/94 (C8) 
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where 



and 



9 9 2- 

2 (j, + to — m F 
sin r = m rr- 5 - 



to. 



- m l( P2 ) ~ m F 



m 



2+ 



(C.9) 
(CIO) 



Since supersymmetry is broken, the complex chiral fields pi do not have holomorphic 
masses. It is necessary to find the combinations of these fields that diagonalize the bosonic 
mass matrix. The eigenvalues of this matrix are 



in 



2 V p 
B 



-rif + m 2 + X 2 + 2/i 2 + py/( v f + m{- \X\ 2 ) 2 + 4p 2 (m 2 + \X\ 2 + 2 Re(X)mi) 



(C.ll) 

where rj = ±1 and p = ±1. The scalar components that diagonalize the boson mass matrix 
are 



B 



C 



-Im(pi+p 2 )cos6' + Im( / 93+ / 94)sin6', 
Im(pi+p 2 ) sin0+Im(/93 + p 4 ) cos/9, £f 
Im(pi - p 2 ) cos 7- Im(p 3 - p 4 ) sin 7 , 

-Im(pi-p 2 )sin7-Im(p3-/> 4 )cos7, 



-Re(/9i — p 2 ) cos 0— Re( / 03 — /04) sin# 
Re(pi — P2) sm9 — Re(ps— pi) cos 
Re(/9i+p 2 ) cos7+Re( / o 3 + / o 4 ) sin 7 

-Re(pi+p 2 ) sin7+Re(p 3 +/?4) cos 7 

(C.12) 



where 



sin 



cos 2 7 



f(m\ 


-X 2 ) + (m 2 f - 


-m 2 +)(m 2 ++- 


m 2 B + + 


m) 


-mf ) 




2(m|+- 


f 2H — \/ 2- 

— j {rrij 


-to 2 +) 






f(m\ 


-X 2 ) + (m 2 f - 


2+\/ 2 — h 

— rrif ){m B - 


- m 2 B - - 


m 2 f 


+ m 2+ ) 




2(m 2 B - 


+ — m 2 B )(m 2 ~ 


-mf) 







The diagonal masses of these fields are 



m 



2+- 



m| 



to 



2+^ 



m 2 c 



to 



B 



m 2 D 



111 



B 



(C.13) 



(C.14) 



We can now evaluate the one loop fermion mass, by using the diagram in figure 10. From 
this diagram one computes the function I(mB,mp), which is 



I(m B ,m F ) = 
and the mass for the fermion is 



16vr 2 



A 2 
lQ g — 



B 



m. 



m 



B 



m, 



■log-T 



(C.15) 



(cos 6 cos t — sin6 , sinr) 2 I(m,^,TO 2 ) + (sin 6 cost + cos0sinr) 2 /(TO j B,TO 2 ) 
(cos 7 cos t + sin 7 sin r) 2 I(mc,TO 2 ) — (sin 7 cost — cos7sinr) 2 /(m,£),TO 2 ) 
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Figure 10: The one loop Feynman diagram associated with the function \(uib, Trip) 



Analogously we can evaluate the mass term acquired at one loop by the scalars Y and Y. 
The mass term in the Lagrangian is 

C = m YY * \Y\ 2 + m\y* \ Y \ 2 + {m Y yYY + cc) (C.16) 

The one loop masses acquired by the pseudomoduli are 

m ry* =m yy* = ((cos 9 cos r — sin 9 sin r) 2 (m 2 +X 2 ) + (cos 9 sin 9+ cos 9 sin T) 2 fi 2 ) K(mA,ni2) 
+ (cos # cost — sin # sin t) 2 J(myi) + (cos6>^ sin6>, sin#— >— cos#, mi-tmg) 
+ ((m^,m B ) ->■ (mc,m D ), 9 -»• 7) 

— ((sin 2 r(m 2 + X 2 ) + /i 2 cos 2 r)K(m 2 7 ~, 7712) + ( cos r — ► sinr, m 2 ^ - — > m 2 ^ )) 

(C.17) 

and 

m y y = 2mA(cos 9 cos r — sin 9 sinr) 2 K(rriA, ^2) + (cos 9 —* sin#, sin 6* — > — cos 9, — > itla — > tub) 
+ ((mi, mg) — > (mc, mo), # —> 7) — sin 2 TK(m 2 F T , m-i) — cos 2 TK(n?p , rri2)) (C.18) 

where the functions K(itibi,iti>B2) and J(ms) are associated with the Feynman diagrams 



of figure |11| and |12[ The computation of these diagrams gives 

. . 1 / A 2 m 2 B1 m? 
K(m jB i,m B 2) =-T^Z?( 1 °g — 2 2 2~~ lo S ~~ 2 



167T 2 \ m^ 2 m B1 — m" m m B2 
.](/»,;) - JL ("a 2 -m^ log ^ (C.19) 



In the R symmetric limit, m-i — * 0, the vev of the scalar pseudomodulus X vanishes. In 
this limit the masses are 

m 2 YY * = m\ Y , = (C.20) 

_ ,u 2 ((l - z^ 2 )((l - e 2 )(^ 2 - 1 - e) log(l - e) + (1 + e 2 )(z^ 2 + e - 1) log(l + e)) - 4eV logv) 

32tt 2 (u 2 -l)(e 2 -{y 2 - l) 2 ) 

where we defined e = -4 and v = At the lowest order in the supersymmetry breaking 
scale this mass reduces to 



m 2 -m 2 - f MM- P (l-^ 2 + 3^-4^1og,) 
myy, - m^, - 327rV ^W - 327r2 ^ 2 (1 _ ^ 2)3 « ( -1) 
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m B2 



K (m B i,m B 2 )= 



m Bi 

Figure 11: The one loop Feynman diagram associated with the function K (771/31,77732) 



\ m B 

J(m B )= - — 

Figure 12: The one loop Feynman diagram associated with the function J(ms) 

where M(v) is a positive function. 
The flavored deformed L 131 model 

In section |] we studied a different dark matter sector in addiction to the superpotential 
(C.l). The superpotential for this DM sector is 



W 3 = Yfafo + Yfofo + /i 2 (</>508 + fo<h) + m 2 <j)7(p8 (C22) 

In this case the squared mass matrices of the fields $5, . . . ,(f)g have to be rotated in a 
diagonal form. This is done by defining the function 

p 2 - X 2 

sin 2 a = (C.23) 

where 



A± = ml + 2p 2 - myjm 2 + 4p 2 ^ ^ 

The diagonal combinations p\ , . . . , pg appearing in W3 are defined by 

05 = — sin aps + cos apg, (pg = sin ap 5 + cos apg 

4>q = — sin ape + cos apj, (pj = cos aps + sin apg (C.25) 
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The fermion mass in this case is 

M^, Y ^ 9 = (cos# cost — sin # sin r) 2 (cos 2 a I(rriA, m/-) + sin 2 a l(rriA, 

+ (sin cos t + cos sin r) 2 (cos 2 a I(ms,m/_) + sin 2 a I(ms,mj + )) 

— (cos 7 cost + sin 7 sin r) 2 (cos 2 a I(mc, + sin 2 a I(mc,m^ + )) 

— (sin 7 cos r — cos 7 sin r) 2 (cos 2 a I(mD,m/_) + sin 2 a I(m.D, 

(C.26) 

where ^ 

m f± = ?»±y/™Z±M (C.27) 
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